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Abstract 

We are concerned with the discretization of a solution of a Forward-Backward 
stochastic differential equation (FBSDE) with a jump process depending on the Brow- 
nian motion. In this part, we study the case of Lipschitz generators, and we refer to 
the second part of this work [11] for the quadratic case. We propose a recursive scheme 
based on a general existence result given in the companion paper [10] and we study 
the error induced by the time discretization. The L^-norm of the error is shown to 
be of the order of the square root of the time step. This generalizes the results for 
Brownian FBSDEs. 
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1 Introduction 



In this paper, we study a discrete-time approximation for the solution of a forward-backward 
stochastic differential equation (FBSDE) with a jump of the form 

Xt = x+ [ b{s,X,)ds+ [ a{s,Xs)dWs+ [ P{s,X,^)dH, , 
Jo Jo Jo 

Yt = g{XT)+[ f{s,X,,Ys,Zs,Us)ds- [ Z^W.-f U^dHs , 
Jt Jt Jt 

where Ht = lr<t and r is a jump time, which can represent a default time in credit risk 
or counterparty risk. Such equations naturally appear in finance, see for example Bielecki 
and Jeanblanc [2], Lim and Quenez [T3], Peng and Xu [15], Ankirchner et al. [1] for an 
application to exponential utility maximization problem and Kharroubi and Lim [10] for the 
hedging problem in a complete market. Our work is divided into two parts. In this part, we 
study the case where the generator / is Lipschitz. The case of a generator / with quadratic 
growth w.r.t. Z is studied in the second part |llj . 

For Lipschitz generators, the discrete-time approximation of FBSDEs with jumps is stud- 
ied by Bouchard and Elie [1] in the case of Poissonian jumps independent of the Brownian 
motion. Their approach is based on a regularity result for the process Z, which is given by 
Malliavin calculus tools. This regularity result for the process Z was first proved by Zhang 
[T6] in a Brownian framework to provide a convergence rate for the discrete-time approxi- 
mation of FBSDEs. The use of Malliavin calculus to prove regularity on Z is possible in [1] 
since the authors suppose that the Brownian motion is independent of the jump measure. 

In our case, we only assume that the random jump time r admits a conditional density 
given W , which is assumed to be absolutely continuous w.r.t. the Lebesgue measure. In par- 
ticular, we do not specify a particular law for r and we do not assume that r is independent 
of W as for the case of a Poisson random measure. 

To the best of our knowledge, no Malliavin calculus theory has been set for such a frame- 
work. Thus, the method used in [1] fails to provide a convergence rate for the approximation 
in this context. 

We therefore follow another approach, which consists in using the decomposition result 
given in the companion paper [10] to write the solution of a FBSDE with a jump as a 
combination of solutions to a recursive system of FBSDEs without jump. We then prove a 
regularity result on the Z components of Brownian BSDEs coming from the decomposition 
of the BSDE with a jump. This regularity result allows to get a rate for the convergence of 
the discrete-time schemes for these BSDEs as in [16] or [1]. 

Finally, we recombine the approximations of the solutions to recursive system of Brow- 
nian FBSDEs to get a discretization of the solution to the FBSDE with a jump. 

We notice that our approach also allows to weaken the assumption on the forward jump 
coefficient. More precisely, we only assume that /3 is Lipschitz continuous, unlike [1] sup- 
posing that (3 is regular and the matrix Id + V (3 is elliptic. 
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As said above, this kind of FBSDEs with a jump appears in finance. The general as- 
sumptions made on the jump time r allow to modelize general phenomena as a firm default 
or simpler as a jump of an asset that can be seen as contagion from the default of another 
firm on the market, see e.g. [5] for some examples. In particular, the approximation of 
these FBSDEs has its own interest, since it provides approximations of optimal gains and 
strategies of the studied investment problems. We study in this part the case of FBSDEs 
with Lipschitz generators, which is related to valuation in complete markets (see [10]) and 
the utility maximization in incomplete markets with compact investment constraints (see 
[11]). The study of the discretization of FBSDEs with a quadratic generator, which are 
related to more general investment problems (see [I]), is postponed to the second part of 
this work. 

We choose to present our results in the case of a single jump and a one- dimensional 
Brownian motion for the sake of simplicity. We notice that they can easily be extended to 
the case of a (i-dimensional Brownian motion and multiple jumps with eventually random 
marks, as in [TU], taking values in a finite space. 

The paper is organized as follows. The next section presents the framework of progressive 
enlargement of a Brownian filtration by a random jump, and the well posedness of FBSDEs 
in this context. In Section 3, we present the discrete-time schemes for the forward and 
backward solutions based on the decomposition given in the previous section. Finally, in 
Sections 4 and 5, we study the convergence rate of these schemes respectively for the forward 
and the backward solutions. 

2 Preliminaries 
2.1 Notation 

Throughout this paper, we let (f2, Q, P) a complete probability space on which is defined a 
standard one dimensional Brownian motion W. We denote F := {J^t)t>o the natural filtration 
of W augmented by all the P-null sets. We also consider on this space a random time r, 
i.e. a nonnegative J-'-measurable random variable, and we denote classically the associated 
jump process by H which is given by 



We denote by D := (I^f)4>o the smallest right-continuous filtration for which r is a stopping 
time. The global information is then defined by the progressive enlargement G := {Qt)t>o of 
the initial filtration where 



£>0 

for all t > 0. This kind of enlargement was introduced by Jacod, Jeulin and Yor in the 80s 
(see e.g. [6], [7] and [5]). We introduce some notations used throughout the paper 



t> . 
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- V(¥) (resp. V{G)) is the a-algebra of F (resp. G)-predictable measurable subsets 
of r2 X M+, i.e. the a-algebra generated by the left-continuous F (resp. G)-adapted 
processes, 

- VAi{¥) (resp. VA4{G)) is the cr-algebra of F (resp. G) -progressively measurable 
subsets of r2 X M_|_. 

We shall make, throughout the sequel, the standing assumption in the progressive enlarge- 
ment of filtrations known as density assumption (see e.g. [H [9l [TO]). 

(DH) There exists a positive and bounded ■p(F) (g) i3(M+)-measurable process 7 such that 

F[Tede\Tt] = lt{0)de , t>0. 

Using Proposition 2.1 in [10] we get that (DH) ensures that the process H admits an 
intensity. 

Proposition 2.1. The process H admits a compensator of the form Xtdt, where the process 
A is defined by 

Xt := , ^lt<r , t>0 . 

F[T>t\j^t] - ' 

We impose the following assumption to the process A. 

(HBI) The process A is bounded. 

We also introduce the martingale invariance assumption known as the (H) -hypothesis. 
(H) Any F-martingale remains a G-martingale. 

We now introduce the following spaces, where a,b E M+ with a < b, and T < 00 is the 
terminal time. 

- S^[a,b] (resp. S^[a,b]) is the set of VAi{G) (resp. PA^(F))-measurable processes 
(yt)te[a,b] essentially bounded 

Il^ll5°°[a6i •= esssup|Ft| < 00. 

te[a,b] 

- S^[a,b] (resp. S^[a,b]), with p > 2, is the set oiPM{G) (resp. PX(F))-measurable 
processes (Vf)fg[a,;,] such that 



sup \Yt\ 

t6[a,b] 



< 00 . 



Hg^[a,b] (resp. H^[a,b]), with p > 2, is the set of V{G) (resp. ■p(F))-measurable 
processes {Zt)t£[a,b] such that 



b . £ 

2 



iZA'^dt 



< 00 . 
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- L'^{X) is the set of P(G)-measurable processes {Ut)telo,T] such that 




) 



llt^ll 



\s\Us\'^ds 



< oo . 



2.2 Forward-Backward SDE with a jump 

Given measurable functions 6 : [0, T] x M ^ M, a : [0, T] x M ^ M, /3 : [0, T] x M ^ M, ^ : M ^ 
M and / : [0, T] x M x R x M x M — M, and an initial condition x G M, we study the discrete- 
time approximation of the solution {X, Y, Z, U) in 5^[0, T] x 5g=[0, T] x i7^[0, T] x L'^{\) to 
the following forward-backward stochastic differential equation 

Xt = x+ [b{s,X,)ds+ [a{s,X,)dW,+ [(]{s,X,-)dHs , 0<t<T, (2.1) 



when the generator of the BSDE is Lipschitz. 

Remark 2.1. In the BSDE (12. 2p . the jump component U of the unknown (Y, Z, U) appears 
in the generator / with the additional multiplicative term 1 — if . This ensures the equation 
to be well posed in iS^[0,T] x H^[0,T] x L'^{\)- Indeed, the component U lives in L'^{\), 
thus its value on (r A T, T] is not defined since the intensity A vanishes on (r A T, T]. We 
therefore introduce the term 1 — if to kill the value of U on (r A T, T] and hence to avoid 
making the equation depending on it. 

We first prove that the decoupled system (12. ip - (12. 21) admits a solution. To this end, we 
introduce several assumptions on the coefficients b, a, (3, g and /. We consider the following 
assumption for the forward coefficients. 

(HF) There exists a constant K such that the functions b, a and /3 satisfy 





< t < T , 



(2.2) 



bit,0)\ + \a{t,0)\ + \/3it,0)\ < K, 



and 



\b{t,x) -b{t,x')\ + \a{t,x) - a{t,x')\ + \P{t,x) - (3{t,x')\ < K\x - x' 



for all (t, X, x') G [0, T] X M X M. 

For the backward coefficients g and /, we impose the following assumption. 
(HBL) There exists a constant K such that the functions g and / satisfy 



\f{t,x,0,0,0)\ + \g{x)\ < K, 
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and 



\f{t,x,y,z,u) - f{t,x,y',z',u')\ < K {\y - y'\ + \z - z'\ + \u - u'\) , 
for all (t, X, y, y\ z, z' , u, u') G [0, T] x R x [R]^ x [R]^ x [M]^. 

In the sequel K denotes a generic constant appearing in (HBL) and (HF) and which may 
vary from line to line. 



Following the decomposition approach initiated by [10], we introduce the recursive system 
of FBSDEs associated with fl23]) -fl2:2D. 

• Find {X\e),Y\e),Z\e)) e 5|[0,T] x S^[e,T] x H^[e,T] such that 

Xl{e) = x + [b{s,Xl{e))ds+ [ a{s,Xl{e))dWs + (3{e,X's-{e))le<t , 0<t<T,(2.3) 
Jo Jo 

= g{x'Ae))+ [f{s,xl{e)X{o),zl{e),o)ds- [zl{e)dw,, e<t<T, (2.4) 

Jt Jt 

for all e G [0,T]. 

• Find {X°,Y^,Z^) G Sl[0,T] x S^[0,T] x H^[0,T] such that 

X° = x + [ [ a{s,X^^)dWs, 0<t<T, (2.5) 

Jo Jo 

= 9{X'^) + J^ f{s,XlY^,ZlY}{s)-Y^)ds- Z'JW, , 0<t<T. (2.6) 

Then, the link between the FBSDE ([211])- dOD and the recursive system of FBSDEs 
(12. 4p and (I2.5p - (I2.6I) is given by the following result. 

Theorem 2.1. Assume that (DH), (HBI), (H), (HF) and (HBL) hold true. Then, the 
FBSDE ([2IID-([22|) admits a unique solution (X, F, Z, U) G 5^[0, T] x 5^[0, T] x iy|[0, T] x 

f Xt = X°li<, + Xi(r)l.<f , 

Ft = YHt<,r + YMT)l^<t , 

* < t I ; ' .2.7) 

= Z°W + Z,i(r)l.<i, 

{ Ut = (F/(t)-F,°)W, 

w/iere {X\e),Y\e), Z\e)) IS the unique solution to the FBSDE ([23D-([23D m S^[0,T] x 
S^[9,T] X i/j^[^,r], for 9 G [0,T], and (X°,F°,Z°) t/ie umgue so/u^^or^ to the FBSDE 
(I23|) ^n 5| [0, T] X 5^ [0, T] x i7| [0, T] . 



Proof. 

Step 1. Solution to under (HF). 
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Under (HF) there exist unique processes X° G S^[0, T] satisfying ([23]), and X\e) e S^[0, T] 
satisfying O for all 9 G [0,T] such that is VM{W) O ;B(M+)-measurable. Then, from 
the definition of H, we check that the process X defined by 



(2.8) 



satisfies (12.11) . We now check that X G Sq[0,T]. We first notice that from (HF), there 
exists a constant K such that 



E 



sup \X^ 

te[o,T] 



0|2 



< K 



Then, from the definition of X^ and X^, we have for all t G [0,T] 

snp\Xl{e)\^ < + |/3(e,X°)|'+ / \b{u, Xl{e))\^ du + snp 

se[e,t] ^ Je se[e,t] 



(2.9) 



a{u,Xl{e))dW^ 



Using (HF) and BDG-inequality, we get 

|2 



E 



sup \xl{e)\ 



< K 1 



E sup \X'^{9)\ 

ue[e,s] 



du 



Applying Gronwall's lemma, we get 



sup ||X^(^)| 

eG[0,T] 



< K 



(2.10) 



Combining ([21]), ([23]) and fl230D . we get that X G >S^[0,r]. Moreover still using (HF) we 
get the uniqueness of a solution to (12. ip in iS^[0,T]. 

Step 2. Solution to ^ under (DH), (HBI), (H) and (HBL). 

To follow the decomposition approach initiated by the authors in [TU] , we need the generator 
to be predictable. To this end, we notice that in the BSDE (12. 2p . we can replace the generator 
(t, y, z, u) H-^ /(t, Xt, y, z, (1 - Et)u) by the predictable map (t, y, z, u) f{t, Xt~ , y, z, (1 - 
Ht~)u). 

Using the decomposition (12. 8p . we are able to write explicitly the decompositions of the 
^7^-measurable random variable giX^) and the P(G) (8)-B(R) ® S(M) (8>S(M)-measurable map 
(cu, t, y, z, u) i-> /(t, Xi- (w), y, z, u(l — H^- (w))) given by Lemma 2.1 in [TU] 

g{XT) = ^7(X°)lT<. + ^?(X^(r))lr>. , 
f{t,Xt-,y,z,{l- Ht-)u) = f{t,y,z,u)tt<r + f\t,y,z,u,T)lt>r, 

with f{t, y, z, u) = fit, X°, y, z, u) and /^(t, y, z, u, 6) = fit, (6*), z, 0), for all (t, z, u, ^ 
[0,T] X R X R X R X M+. 

Suppose now that (DH), (HBI), (H) and (HBL) hold true. Then, from Theorem C.l 
in [in], the BSDE ([23]) admits a P(F) ® -B([0, T])-measurable solution (F\ Z^) and the 
BSDE (12. 6p admits a solution iY^,Z^). Using Proposition 2.1 in [13], we obtain 

\\YHe)\\.^,a^. + WzH 



< K 



for all e e[0,T], and 



ll'^°ll5°°[0,T] + ll^°ll//2[0,T] — ^ ■ 

We can then apply Theorem 3.1 in [10] and we get the existence of a solution to (12. 2p in 
S^[0,T]xHl[0,T]xL\X). 

Let {Y,Z,U) and {Y',Z',U') be two solutions to (Q in S^[0,T] x H^[0,T] x ^^(A). 
Since f(t, x, y, z, {1 — Ht)u) = f{t, x, y, z, 0) for all t G (r AT, T] and A vanishes on (r AT, T], 
we can assume w.l.o.g. that Ut = Uj. = for t G (r AT, T]. Then, from (DH), (HBI), (H) 
and (HBL), we can apply Theorem 4.1 in [10] and we get that Y < Y'. Since Y and Y' 
play the same role, we obtain Y = Y'. Identifying the pure jump parts of Y and Y' gives 
U = U'. Finally, identifying the unbounded variation gives Z = Z'. □ 

Throughout the sequel, we give an approximation of the solution to the FBSDE (12. ip - 
(12. 2p by studying the approximation of the solutions to the recursive system of FBSDEs 
(E^D-dSaD and CT-d^. 

3 Discrete-time scheme for the FBSDE 

In this section, we introduce a discrete-time approximation of the solution {X, Y, Z, U) to 
the FBSDE (I2.ip - (l2.2p based on its decomposition given by Theorem 12.11 

Throughout the sequel, we consider a discretization grid n := {to, . . . of [0,T] with 
= to < ti < . . . < tn = T . For t G [0, T], we denote by 7i{t) the largest element of vr smaller 
than t 

ii(t) := max {ti , i = 0, . . . , n | < t } . 

We also denote by |7r| the mesh of vr 

|7r| := max | tj+i — ti , ^ = 0, — l}, 

that we suppose satisfying |7r| < 1, and by AWl' (resp. At^) the increment of W (resp. the 
difference) between ti and tj^i: AW[ := Wt^ - Wti_^ (resp. At^ := ti - for 1 < z < n. 

3.1 Discrete-time scheme for X 

We introduce an approximation of the process X based on the discretization of the processes 
and X^. 

• Euler scheme for We consider the scheme X^''^ defined by 

*° ~ (3 1) 

Xj- = Xj^+6(t,„i,Xj:^)Atr + a(t,_i,X°:-JAW7, l<i<n. 
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• Euler scheme for X^. Since the process X^ depends on two parameters t and 9, we 
introduce a discretization of X^ in these two variables. We then consider the following 
scheme 

^i''(7r(6')) = X + /3(to,a;)l^(e)=o , 

4'-(vr(^)) = xl::^{Tx{e)) + h{u.,,xl;:Me)))/\t^ + a{u.,,xl^^^^^ (3.2) 

+ /3(t,_i,X;':^(7r(e)))lt,=,(,) , l<i<n, 0<e<T. 

We are now able to provide an approximation of the process X solution to the FSDE (12.11) . 
We consider the scheme X'^ defined by 

X: = Xl'^^t,^, + Xl';^inir))t,^, , 0<t<T. (3.3) 

We shall denote by {J^i''^}o<i<n (resp. {J^i'^{d)}o<i<n) the discrete-time filtration associated 
with (resp. X^'^) 

:= CT{X?;\J<^) 
(resp. J-/'-(^) := a {Xlf (9) , j < z)) . 

3.2 Discrete-time scheme for (Y, Z, U) 

We introduce an approximation of (y, Z) based on the discretization of (F*^, Z°) and (y^, Z^). 
To this end we introduce the backward implicit schemes on vr associated with the BSDEs 
(12. 4p and (12. 6p . Since the system is recursively coupled, we first introduce the scheme 
associated with (12. 4p . We then use it to define the scheme associated with (12.61) . 

• Backward Euler scheme for (F^, Z^). We consider the implicit scheme {Y^''^ , Z^''^) defined 
by 

f Y^^-{n{e)) = giX'^^m)) , 
yt\im) = E,^lf)K-(7r(^))] +/(t,_i,X,^-(vr(0)),ri';(^ 

(3.4) 

where E]'' = E[ . \TI'^{s)] for < z < n and s G [0,T]. 

• Backward Euler scheme for {Y^,Z^). Since the generator of (12. 6p involves the process 
(l^^(t))^gjQ -Pj, we consider a discretization based on Y^''^ . We therefore consider the scheme 
JyCTT^ ZO''^) defined by 



yt\ = Et,K-]+r(t,_i,x°-,it:,^S::^)Air 



(3.5) 
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where E° = E[ . | J-"^ ''^] for < z < ra, and is defined by 

l,7r 
r(t) 



for all {t, X, y,z) e[0,T]xRxRx R. 

We then consider the following scheme for the solution (Y, Z, U) of the BSDE fl2.2p 

Yt"" = i;°5W + F\-(7r(r))l,>., 

Z: = Z;'J1,<, + (7r(r))l,>. , (3.6) 
Ut = K(J(vr(t))-F:'-)W, 

for t e [0,T]. 

4 Convergence of the scheme for the FSDE 

We introduce the following assumption, which will be used to control the error between X 
and X"". 

(HFD) There exists a constant K such that the functions b, a and /3 satisfy 

\h{t,x) -b{t\x)\ + \a{t,x) - a{t\x)\ < K\t-t'\^^, 
\p{t,x)- (5{t\x)\ + \a{t,x)-a{t',x)\ < K\t-t'\, 

for all (t,t',x) G [0,T] x [0,T] x R. 

In the following we provide an error estimate of the approximation schemes for X° and 
X^ which are used to control the error between X and X'^ . 

4.1 Error estimates for and 

Under (HF) and (HFD), the upper bound of the error between X^ and its Euler scheme 
X^''^ is well understood, see e.g. [12], and we have 



E 



sup |X°-X;j |' < J^lvrl, (4.1) 



for some constant K which does not depend on tt. 

The next result provides an upper bound for the error between X^ and its Euler scheme 
X^'^ defined by flX^ . 



Theorem 4.1. Under (HF) and (HFD), we have the following estimate 

< K\'k\ , 



sup E 

9g[0,T] 



sup \xl{e)-xl';^{7^{e)f 



for a constant K which does not depend on tt. 
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Proof. Fix 9 G [0,T], we then have 



E 



sup \xl{e)-xl';{'n{e))\ 



< 2 E 
+ 2 E 



sup \xl{e) - xl{^{e))\ 



sup \Xlinie)) - Xl';^ini9))\' . (4.2) 



t€[e,T] 



We study separately the two terms of the right hand side. 

Since 71(6) < 6 < t, we have by definition Xj(7r(^)) = for all s G [0,7r(^)), and = 
X^ for all s G [0, 6), which implies 

Xl{e)-Xlinie)) = f b{s,X'^)ds+ f a{s,X'^)dWs + ^{e,X'e)+ fK',Xli9))ds 

Jirie) Jn(e) Je 



<y[s,xl{e))dw,- p{T,{e),x%^) - / b{s,xl{T:{e)))ds 



a 



{s,xl{ii{e)))dWs 



for all t G [e,T]. 

Hence, there exists a constant K such that 



\Xl{e)-Xl{7,{e))\^ < k{ f b{s,X^)ds\ f b{s,Xl{n{9)))d, 

Jwie) Jn(e) 

b{s,Xl{e))-b{s,Xl{nie))) ds 
a{s,X^)dWs 



AG) 

t 



a 



[s,xl{Tx{e)))dWs 



a{s,Xl{e))-a{s,Xl{7,{e))))dWs 



+ |/3(^,X°)-/3(vr(^),X°(,))|' 



(4.3) 



From (HF) and (HFD), we have 



E|/3(e,X°)-/3(vr(^),X°(,))f < A'(|7rp + E|X°-X°(,)|^) 



x^O 



We have from (HF) and ^B) 



E 



a{s,X^^)dWs < K\n\ 



Trie) 



which implies in particular E|Xg — X^^^-jp < K\7t\ and hence 

E\/3ie,X^,)-/3inie),X%^)\'<K\7r\. 
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We have also from (HF) and (^M) 



E 



7r(0) 



6(.,Xl(7r(^)))d' 



a{s,Xl{Tr{e)))dWs < K\7t 



2i 



Combining these inequalities with (14. 3p . (HF) and BDG- inequality, we get 



E 



sup \xl{e) - xl{n{e))\ 



< K E 



sup \xl{e) - xl{7r{e))\ 

ueie,s] 



ds 



7r 



Applying Gronwall's lemma, we get 

E 



sup \Xl{e)-Xl{Tr{6))\ 



< A'IttI . 



(4.4) 



To find an upper bound for the term E[supjg[g |X^-'^(7r(6')) — X^'^^^{it{6))\'^] we introduce the 
scheme X'^{7r{6)) defined by 



Xl{Ae)) = Xl_^{n{e)) + b{U_,,Xl_^{n{e)))At^ + a{U_,,X^^^^^ u>n{e). 
We have the inequality 

e[ sup \xl{7r{e))-xl^;{ix{e))\'] < 2e 



+ 2 E 



Since X^{7r{e)) is the Euler scheme of X\7r{e)) on [7r(^),r], 
(HFD) (see e.g. [12]) 



■ sup \xl{^{9)) - x:^,^{n{e))\' 
' sup \x:^,^i7rie))-xl;{n{e)) 

we have under (HF) and 



(4.5) 



E 



< k(^i + e 



\Kie)im)\ 



In] 



for 



sup \xl{n{e)) - x:^,^{7r{e))\' 

some constant K which neither depends on tt nor on 6. From (12.101) . we get 

e\ snp \Xl{nie)) - X:^,){7rie))\'] < A>| , 



for all e G [0,T]. 
We now stuc 



(4.6) 



n{t)i^i^)) - Xj)(^(^))r]- We first notice that 



e[ sup |x;(,)(7r(0))-X]f)(vr(^))|'j = e[ sup |x;(,)(7r(^)) - X]f)(vr(^))f" 
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Hence we can work with the second term. From the definition of X'" and X ''^ , we get 



sup |X;(„)(7r(e))-Xij(7r(^))|^< 

MG[7r(e),t] 



K2 




1 + 


L(e) 



Ki\Xl^oM9))-Xll-j7ri9))\'+ / &(7r(s), X;(,)(vr(^))) - 6(vr(.), X^- (7r(e))) 



ds 



+ sup 

ue[7r{e),t] 



■k{u) 



cr 



(7r(.),X;(,)(vr(e)))-a(7r(s),Xij)(7r(^))))dW^. 



Then, using (HF) and BDG-inequahty, we get 



E 



sup _ \x:^uMO)) - Xl'^^^inie))]'] <K(E\Xl^,^inie)) ~ Xl';^{iTi9)) 



«e[7r(6»),t] 



+ / E 

ln{e) 



sup X;(„)(7r(^))-X;^- (7r(^)) 



ds 



From Lipschitz property of (3, we have 



E|x:(,)(vr(e)) -X]5)(vr(^))|^ = E|X°(,) + /3(7r(^), - X^J- - /3(vr(e), 



< ke|x°(,)-x:;(,^| . 



This last inequahty with (14. ip gives 

E|x:(,)(vr(^))-X^5)(vr(^))p < 



Applying Gronwall's lemma, we get 
E 



sup X;(,)(7r(^))-X^J(7r(e)) 
^te[7r{e),T] 



< K\n\ 

Combining (H^j), (jO]), (H^]), dSS]) and (jlZD, we get the result. 



(4.7) 



□ 



4.2 Error estimate for the FSDE with a jump 

We are now able to provide an estimate of the error approximation of the process X by its 
scheme X'^ defined by (13.31) . 

Theorem 4.2. Under (HF) and (HFD), we have the following estimate 

|2" 



E 



sup \Xt-X^\ 

te[o,T] 



< K IT 



for a constant K which does not depend on vr. 
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Proof. From the definition of X'^, (DH) and (14.11) we have 



E 



sup \Xt-X;^\ 
te[o,T] 



< E 



< E 

< K 

From Theorem 14.1^ we get 



sup 






sup 

^te[o,T] 







+ E 


sup 1 




^te[T,T] 







+ [e 


sup 


Jo 





d9 



TT 



sup E 

ee[o,T] 



sup \xl{e) - xX^{n{e))\ 



E 



sup \Xt-X^\ 

t&[0,T] 



< KM . 



□ 



5 Convergence of the scheme for the BSDE 

To provide error estimates for the Euler scheme of the BSDE, we need an additional regu- 
larity property for the coefficients g and /. We then introduce the following assumption. 

(HBLD) There exists a constant K such that the functions g and / satisfy 

\g{x) -g{x')\ + \ f{t,x,y,z,u) - f{t' ,x' ,y, z,u)\ < K{\x-x'\ + \t-t'\^) , 
for all {t, t\ X, x', y, z, u) G [0, Tf x [M]^ x R x M x M. 

We are now ready to provide error estimates of the approximation schemes for and 
and then for 



5.1 Regularity results 

In this part, we give some results on the regularity of the processes and We denote 
7"° := , < s < t} and := a{X\{e) , ^ < s < t}. 

Proposition 5.1. Under (HF), (HFD), (HBL) and (HBLD), there exists a constant 
K such that 



E 



\zl{e)-zl^,^{efdt 



< K l + E 



\X'ei0)\ 



TT 



(5.1) 



for all 6 E TT. 



Proof. We first suppose that b, a, f and g are in C^. Let us define the processes A and M 
by 



A, := exp( f dyf{Ql{e))dr) , 
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and 



Mrdj{elie))dWr 



where 01(9) := {r, XX9),Y^\9), Zl{9),0). We give classically the link between V^X/(^)(:= 
dXl{9)/dX^{9)) and {DsXf{9))e<s<t the Malliavin derivative of Xl{9). Recall that X\9) 
satisfies 

Xl{9) = X^g{9) + [b{r, Xl{9))dr + [a{r, Xl{9))dWr , 9<t<T . 
Je Je 

Therefore, we get 

V'Xl{9) = 1+1 d^h{r,Xl{9))V'Xl[9)dr+ ! d^a{r,Xl{9))V' Xl{9)dWr , 9<t<T 
J 9 Je 

and for 6* < s < t 

DsXl{9) = a{s,Xl{9))+ [ d^h{r,Xl{9))DsXl{9)dr+ j d^a{r,Xl{9))DsXl{9)dWr . 

J s J s 



Thus, we have 



D,Xl{9) = V'Xl{9)[V'Xl{9)]-'a{s,Xl{9)) . 



(5.2) 



Using Malliavin calculus we obtain that a version of Z^{9) is given by {DtYl{9))^^yQrj,y By 
Ito's formula, we get 



KtMtZ]{9) = ¥.\MT[kTVg{Xl{9))DtXl{9) + dj{Ql{9))DtXl{9)Kdr) J'l{9) 
for t e [9,T]. Using (Q, we get 

KMtZl{9) = E^MT(^AT'^g{X^{9))V^X^{9) + FrArdr^ J^l{9)j[V^Xl{9)]~'(T{t,Xl{9)) 
with Fr := dj{el{9))V^X^{9). This implies that 

AtMtZl{9) = (E[G\Tl{9)]-Mt f FrArdr)[V'Xl{9)]-'a{t,Xl{9)) , 



with G := MT(^ATVg{X^{9))V^X^{9) + FrArdr^ Since 6, a, / and g have bounded 
derivatives, we have 

E[|G|P] < oo, p>2. (5.3) 
Define := E[G| J^^^(6')] for r G [9,T]. From (15.31) and Doob's inequality, we have 

||"^||5p[e,T] < oo , p > 2 . (5.4) 
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Hence, there exists a process </) such that 

= E[G\J^^{e)]+ [ <PudW^, re[^,T], 
Je 

and 

||0||//P[0,T] < OO , p > 2 . 

We define Z by 

J 6 

By Ito's formula, we can write 

Zt{e) = Ze{9)+ [ al{e)dr+ [ al{e)dWr , e<r<T. 
Je Je 

Since b, a, f and g have bounded derivatives, we get from (15 ■4p 



and 



sup ||^(^)||5prerl < ) P > 2 

ee[o,T] 



sup (\\a^i6)\\HPie,T] + \\a^iO)\\HP[e,T]) < oo , p>2. 



2 ) 



6»G[0,r] 

We now write for t G [ti, tj+i 
with 

ily.= E[\z,ie)\'\a{t,xlie)) - a{u,xl{e))\'] . 

We give an upper bound for each term. 



< KE 



E[\z,ie) - z,xo)\'\Tlie)]\a{u,xlie))\ 

'^^\\al{e)\^ + \a'',{e)\^)dr sup \a{t,Xl{e))f 



which imphes 



Il^dt < K\tt\E 



t€[e,T] 



'^^\HW + HW)dr sup \a{t,Xl{e))\ 



therefore we have 



n-l 

j=o, ti>e'^^i 



ti+i 



1} At < K\tt\E 



{\al{e)\^ + \al{e)\^)dr sup \a{t,Xl{e))\ 



t€[e,T] 
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From Holder's inequality and (HFD) and (HE), we have 



i=o, ti>e 
Using ( I5.6p . we get 



,2|4 



{9))dr 



1 +E 



sup |X/(^)| 

t&[e,T] 



n-l 



i=o, ti>e 



ti + l 



Il.dt < K\tx\[1 + ¥. sup 



< K\tt\ 1 + E 



We get from ([53]), (HFD) and (HE) 



It, < KiE 



\z,-z,X\xlxe)\ 

Arguing as above, we obtain 

"ti + l 



E 



\xl{e)Zt-xl{e)Zu\ 



TT 



n-l 



i=0, ii>0 " 



E 



rft < AIttIE 



sup {l + \Xlm') 



e<t<T 



Moreover, from Ito's formula, X^{9)Z is a semimartingale of the form 

Xl{d)Zt = Xl{e)Ze+ [ aldr+ [ a^dWr , 



where | |_h-2[6»,t] + I l^^l I_ff2[5i^r] < K{l + E[\X^{e)\^]4). Therefore, we have 



E 



\xl{e)z,-x}{e)z,^ 



< KE 



;?1|2 , |~2|2 



\a. 



'')dr 



'-Jt 



which implies 



(5.7) 



(5.8) 



j ^^E[\xl{e)z,-xl{e)Zu\'\ < k\t,\e 

i=o, ti>e 

Using (15. 7p and (15. 8p we get the result. 

When 6, a, /3, / and g are not in C^, we can also prove the result by regularization. We 
first suppose that / and g are in C^. We consider a density q which is on M with a 
compact support, and we define an approximation (6^, cr^, /3^) of (6, a, /3) in Cl by 



(6SaS/3^)(t,x) = - / (6,a,/3)(t,x')g 



X — X 



rfx' , (t, x) G [0, T] X M . 



We then use the convergence of {X^'\9) ,Y^'%e) , Z^'^O)) to {X^{e),Y^{e), Z^{e)) and we 
get the result. Next we assume that / and g are not Cl and we consider for that and g"^ 
which are defined as previously and we get the result. □ 

Using the link between X^ and (6'), we obtain that the bound (15.11) is actually uniform 



in e. 
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Corollary 5.1. Under (HF), (HFD), (HBL) and (HBLD), there exists a constant K 
such that 



E 



\zl{e)~zl^^{e)\'dt 



< K\'r\ 



(5.9) 



for all 6 ^ IT. 



Proof. Since is a Brownian diffusion, we have for any p > 2, from (HFD) and (HF), 
tliat 



E 



sup 

te[o,T] 



\x?\p 



< oo 



We notice tliat from the Lipschitz property of /3 we have 



E 



E 



< K 1 + E 



sup 

te[o,T] 



< oo 



Combining this result with (15.11) . we get (15.9 
We now study the regularity of Z^. 



□ 



Proposition 5.2. Under (HF), (HFD), (HBL) and (HBLD), there exists a constant 
K such that we have 



E 



" ^n(t) I dt 



{t)\ 



< KM . 



Proof. The proof is similar to the previous one. The only difference is that the BSDE (12. 6p 
involves Y\ We denote 6° = {r , , {r) - Y^^). We first suppose that 6, a, /3, / 

and g are in Cl- We recall that 



Z'idWs 



Therefore, for 0<r<t<T, we have 



DrY,' = VgiX^)DrX^ + J^^ [d.j{e',)DrX',+{dy-d^)f{e',)DX 



where DX^, DY^, DZ^ and DY^r) denote the Malliavin derivatives of Y^^, Z° and 
Y^{r) for r G [0,T]. Using Malliavin calculus, we obtain that a version of Z^ is given by 
{DtYl^)t(z[o^T]- By Ito's formula, we get 



AtMtZt 



E 



Mt(AtV(7(X°)AX°+ / (9,/(0°)AX° + 9,/(e°)Ai;'W)A.dr 



■J t 



where At := exp(/p (9^^ - (9„)/(G°)(ir) and Mt := 1 + /g M,9,/(e°)dPK. Denote by VXf : = 
g| and VX^{e) := for < t < ^ < T. We then have for r < s < T 

= (1 + 9./3(s,X°))A.X° = (1 + 9,./3(s,X°))VX°a(r,X°)[VX°]-i , 



thus we can see that DrXl{s) = VXj(s)cr(r, X°)[VX°] ^ Therefore, we get by writing the 
SDEs satisfied by {DrXl{e))se[e,T] for r<e, and (VXi(^)),g[0,T] 

DrXHe) = vxHe) [VX°] ~V(r, X°) , r < ^ < s . 

Writing the BSDEs satisfied by {DrY}{e))s(.[e,T] ioi r < 9 and {WY}{9))se[e,T], and using the 
previous equahty, we get 

DrY,\s) = VYXs)[VX^y'a{r,X^) , s<e. 

This implies 

with Fr := 9^./(e°)VX0 + 9„/(e°)Vi;i(r). We can write 

AtMtZt = (e[G|J-°]-^ MiF,A,dr)[VXO]"V(t,XO) , 

with G := M'r(A2-V5'(X^)VX^ + F^A^dr). Since 6, a, / and g have bounded derivatives, 
we have 

E[|G|P] < oo, p>2. (5.10) 
Define := E[G'|J-'°] for r G [0,T]. From flS.lOp and Doob's inequality, we have 

||"^IUf[o,r] < oo , p>2. (5.11) 
Hence, there exists a process such that 

- E[G']+ / ct>.udW^, r G [0,T] , 

^0 



rrir 



and 

||0||hp[o,t] < oo , p>2 . 

We define Z by 

Zt := (A^Mi)-^ (mi - Mt ^ F,A,rfr) [VX°] t G [0, T] . 
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By Ito's formula, we can write 

Zt = Zo+ [ aids + [ aldWr , t G [0, T] . 
Using the fact that 6, a, / and g have bounded derivatives and f lS.lip . we get 



and 



||a^||Hp[o,T] + ||«^||hp[o,t] < oo , p>2. 
We now write for t G [tj, tj+i) 

with 

Ily.= E[|Z,-4|V(t.,^°)r] , 
J,^,:= EOZ,|V(t,X°)-a(t.,X°)p]. 
As previously we give an upper bound for each term. 



lit < 



{\al\'^ + \a'^^\'^)dr sup 
te[o,T] 



From Holder's inequality and Lipschitz property of a, we have 



n-l 



1=0 



^ / Il.dt < K\n\E {\al\' + \a'X)dr 



JO 



1 + E 



sup |X, 
^te[o,T] 



0|4 



Using fl5.12l) . we get 



i+i 



i=0 



From (HFD) and (HF), we get 



< 



k(e 



E 



\x',z,-xlzA 



Arguing as above, we obtain 



n-l 

E 

Moreover, X^Z is a semimartingale of the form 



TT 



/ E 


'\Zt-z,X\xl\\ 


dt < /i|7r|(^l + E 


sup 




Iti 




Ue[o,r] J 





X^Zt = X^Zo+ / &ldr+ / a^diy, 
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where | |//2[qj^] + 1 |_f/2fQj^j < K and we have 



[0,T] 



E 



+ \al\^)dr 



which imphes 



E 

■t=0 



U+1 



E 



< K\ti\ 



When 6, a, / and g are not C^, we can also prove the result by regularization as for 
Proposition 15.11 □ 



5.2 Error estimates for the recursive system of BSDEs 

We first state an estimate of the approximation error for (F^, Z^). 

Proposition 5.3. Under (HF), (HFD), (HBL) and (HBLD), we have the following 
estimate 



sup < sup E 

6»G[0,T] te[0,T] 



E 



zm - z];;,{n{e))\\s 



< K\ti\ 



for some constant K which does not depend on vr. 
Proof. Fix e G [0,T] and t G [e,T]. We then have 



E 



< 2 E 



\Y,\e)-Y,\m)\ 



+ 2 E 



(5.13) 



We study separately the two terms of right hand side. 

Define 5Xl{e) := Xm - X} {i, {9)) , 5Y,\e) := Y,\9) - Y,\7i {9)) and 6Zl{9) := Z\{9) 
Zj-{n{9)). Applying Ito's formula, we get 



\SY^ 



2 / 5YX9) f{elini9)))-f{Qli9)) 



ds 



+ 2 [ 5Y}{9)5Zl{9)dWs+ I \6Zl{9)\^ds 
Jt Jt 

where el{9) := (s, Xj(^), Zj(^), 0). From (HBL) and (HBLD), we get 



E[\6Y,\9)\'] < K(E[\6Xm\']+^ / m\9)\\6Xl{9)\ds 



T 



E 



T 



\5Yl{9)\'ds 



E 



\5Y}{9)\\5Zl{9)\ds 



E 



T 



\5Zl{9)\'ds 
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Using the inequality 2ab < a'^/rj + rjb'^ for a, 6 G M and 77 > 0, we can see that 



E[\6Y,\e)\']+E r\5Zlie)\'ds] < K(E[\6Xm\']+ rm5Y,\e)\']ds 



+ E 



Jt 



From fl4.4p and Gronwall's lemma, we get 



(5.14) 



(5.15) 



We now study the second term of the right hand side of fl5.13p . Using the same argument 
as in the proof of Theorem 3.1 in [5], we get from the regularity of given by Corollary 



E 



\Y,\n{e))-Yli;^{n{e)) 



< K\'k\ . 



(5.16) 



This last inequality with (15.131) and fl5.15p gives 

sup { sup E\\Y,\e)-Ylg{ix{e))\ 



< JilTrl . 



We now turn to the error on the term Z^{6). We first use the inequality 



E 



T 



\zm-zl';^inie))\'dt 



< 2 E 



T 



\zlinie))-zl;^inie))\'dt 



+ 2 E 



T 



\6zl{e)\'^dt 



(5.17) 



Using flsni) and flSTTSj) with t = 6*, we get 



E 



\6zl{e)\^ds 



< K\7r\ . 



(5.18) 



The other term in the right hand side of (15.170 is the classical error in an approximation of 
BSDE. Therefore, using Corollary 15.11 and (I5.16p . we have 



E 



T 



\Zl{7r{9))-Zl';^ini9))\'dt 



< K\ 



TT 



(5.19) 



Combining flSTTj) . fl518D and fl5A9D . we get 



E 



\zKe)-zl';j7r{e))\'dt 



< K\it\ 



□ 



We now turn to the estimation of the error between [Y^, Z^) and its Euler scheme (13.50 . 
Since this scheme involves the approximation Y^''" of Y^, we first need to introduce an 
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intermediary scheme involving the "true" value of the process Y^. We therefore consider the 
scheme [Y^''^, Z^''^) defined by 



(5.20) 



70,7r 



At. 



n . 



Using the regularity result of Proposition 15.21 and the same arguments as in the proof of 
Theorem 3.1 in [3], we get under (HF), (HFD), (HBL) and (HELD) 

-T 



sup E 

*6[0,T] 



E 



I - I dt 



< Klnl 



(5.21) 



With this inequality, we get the following estimate for the error between {Y^, Z^) and the 
Euler scheme (13. 5p . 

Proposition 5.4. Under (HF), (HFD), (HBL) and (HBLD), we have the following 

estimate 

f-T 



sup E 

*e[o,T] 



1^0 _ VA0,7r|2 



E 



I - Zl'^t) I d^ 



< K\ 



TT 



for some constant K which does not depend on n. 
Proof. We first remark that 



sup E 

t6[0,T] 
T 



E 



■-JO 



< 2 sup E 

tG[0,T] 

< 2E' 



+ 2 sup E 

tG[0,T] 

2E' 



■-JO 



/70,7r '70,77 I 2 7, 
7r{t) 7r(t) I 



Using (|5.2ip, we only need to study sup^g^.T] -F^^'^'jl^] ^^'^ ^Uo -^!(t)l^^^]- To 

this end, we need to introduce continuous schemes for all < i < n — 1. Since EdF^^'^^p] < oo 
and E[|y;°'^|2] < oo for all 1 < i < n, we deduce, from the martingale representation theorem, 



that there exist square integrable processes Z^'"^ and Z^' such that 



~0,7r 



We then define 



Y 



0,7r 



Y 



0,7r 



Z^fdW, 



z'fdw. 



rr = 









rOjTT -1^0,71 r^O,TT 



~0,7r 
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for t G [ti, tj+i). Let i E {0, . . . ,n — l} he fixed, and set SYt := Y^'^ — Y^'^ , SZi := — 

5|^:= and := f{u, xJ^ i^;'-(t,)-i^°-)-/(t„ xJ^ y;°'^ 

Y'j.''^) for t G By Ito's formula, we compute that 



i+i 



•i+i 



E\6Z^\'ds-E\6Yt^\' = 2 / E[6YJf.,]ds , < t < . 



Let a > be a constant to be chosen later on. From the Lipschitz property of / and the 
inequality 2ab < aa? + b'^ /a, we get 



At < a 



ti + l 



E\5Ys\^ds + 



a 



E 



\5Yu\' + \5Z,\' + \Yl{U)-Y,Y{U)\' 



ds 



Using Proposition 15.31 we get 



At < a E\SYs\'^ds + —\iT\E\SYt.\'^ + — E\5Zi\'^ds + — 
Jt a " a Jt a 



We can write 



E|511|2 < E\SYt^^f + j^^''\\6Zfds < a j^^^' E\5Ys\'^ds + Bi, (5.22) 



where 



K 



K. 



K. 



n^Yt^+.V + — kl E\5Z,Y + -kl E\5Yty + — Ivr 



a 



a 



a 



By Gronwall's lemma, this shows that E|(51^p < fijc"''^' for ti < t < tj+i, which plugged in 
the second inequality of f l5.22p provides 



E\6Yt\' 



j^^^' E\6Zfds < Bi(^l + a|7r|e"l''l 



(5.23) 



Interpreting Z^^"" (resp. Z^^"") as the projection of Z^''^ (resp. Z^' ) in on the set 

of constant processes, we have 



u 



E\SZi\^ds < / E\SZfds . 



(5.24) 



Applying (15.231) for t = and a = 2K, and using the previous inequality, we get 

E\6Ytf + h{7i) r^\\5Zfds < k2{7r)E\6Yt^^f + k3{TT)\TTW 0<i<n-l 
Jti 



where kiin) 



_i-Kl7rle2^W_ 

l_M_X|7r|2e2^l'^l ' 



= , and k 



l_M_K|7r|2e 



^3 TT 



Since 



l_M_/^|^|2e2/fh 



for small |7r| we have ki{'iT) > 0, we get 



E\6Ytf < k2{n)E\6Yt^J^ + k3{n)\7iW 0<i<n-l, 



24 



for |7r| small enough. 

Iterating this inequality, we get 



Since ^2(71") > 1 and 5Yt^ = 0, we get for small |7r| 

E|5Ft,P < \TT\k3{7r)k2{7r)^ < K\7r\ , 0<i<n 

which gives 

< K\Tr\ . 



(5.25) 



sup E 

t&[0,T] 



Summing up the inequality fl5.23p with t = ti and a = 2K and using fl5.24p . we get 

n— 1 



1 



K|7r|e 



-I)/ E|Z;j)-Z;jj2ds < 2A>|e2^W^E|5Ffj2 + (l + 2ir|7r|)E|5Ff„ 

^ n— 1 

+ (- + i^lvrle^^Wj (^1^1 + |vr| ^E|5Fi, 



Using f l5.25p . we get for |7r| small enough 
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5.3 Error estimate for the BSDE with a jump 

We now give an error estimate of the approximation scheme for the BSDE with a jump. 

Theorem 5.1. Under (HF), (HFD), (HBL) and (HBLD), we have the following error 
estimate for the approximation scheme 

sup E 

tG[o,r] 

for some constant K which does not depend on it. 
Proof. 

Step 1. Error for the variable Y. Fix t G [0,T]. From Theorem 12. II and (13. 6p . we have 
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Using (DH), we get 
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Using Propositions 15.31 and I5.4[ and since t is arbitrary chosen in [0,T], we get 
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Step 2. Error estimate for the variable Z. From Theorem 12. II and (13. 6p . we have 
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Using (DH), we get 
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From Propositions 15.31 and 15. 4[ we get 
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Step 3. Error estimate for the variable U. From Theorem 12. II and (13.61) . we have 
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Using (HBI), we get 
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Combining this last inequahty with Propositions 15.31 and 15. 4[ we get the result. 



□ 



Remark 5.1. Our decomposition approach allows us to suppose that the jump coefficient 
/3 is only Lipschitz continuous. We do not need to impose any regularity condition on /3 
and any ellipticity assumption on Id + V/3 as done in [4] in the case of Poissonian jumps 
independent of W . 
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